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Abstract 


Paraboloid antenna surfaces suffer performance degradation due to structural deforma- 
tion. A first step in the prediction of the performance degradation is to find the best-fit 
paraboloid to the deformed surface. This paper examines the question of whether rigid 
body translations perpendicular to the axis of the paraboloid should be included in the 
search for the best-fit paraboloid. It is shown that if these translations are included the 
problem is ill-conditioned, and small structural deformation can result in large transla- 
tions of the best-fit paraboloid with respect to the original surface. The magnitude of 
these translations then requires nonlinear analysis for finding the best-fit paraboloid. On 
the other hand, if these translations are excluded, or if they are limited in magnitude, the 
errors with respect to the restricted “not-so-best-fit” paraboloid can be much greater than 
the errors with respect to the true best-fit paraboloid. 


Introduction 

Paraboloid antenna surfaces suffer performance degradation due to structural deforma- 
tion. A first step in the prediction of the performance degradation is to find the best-fit 
paraboloid to the deformed surface. The process of finding this best-fit paraboloid has 
received some attention in the past, (e.g. Refs. 1,2) but there is no clear agreement on 
a procedure that should be followed. In particular, questions that arise are whether it is 
acceptable to change the focus length in choosing the best-fit paraboloid and which rigid 
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body motions should be included in moving from the original paraboloid to the best-fit 
one. The present paper attempts to shed some light on this second question. 

The choice of rigid body modes to be considered is associated with ill-conditioning of 
the numerical process of obtaining the best-fit paraboloid. If Z denotes the paraboloid 
axis symmetry, then the ill-conditioning is associated with translations in the X and y 
directions. Because antenna surfaces are typically shallow paraboloids, finding X and y 
translations required to move the original paraboloid to the best-fit one leads to an ill- 
conditioned set of equations. It is possible to eliminate these translations by, for example, 
setting them to be equal to the corresponding translations at the apex. However, it is not 
clear how much is lost in terms of the root mean square (rms) surface error. This paper 
shows that not including these translations can indeed result in substantial increase in 
rms errors, but that to include them one must resort to complicated and costly nonlinear 
calculations. This is demonstrated first by the simpler case of a best-fit parabola. 


Best-Fit Parabola 


The undeformed shape of the parabola is given as 

yi = 


The distortions in the X\ and y\ directions are given by £ and T) (see Figure 1) so that 
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a point A moves to position A'. The best-fit parabola is given as 

V = cx 2 (2a) 

c = a + b (2b) 

The two coordinate systems shown in Figure 1 have unit vectors i\ and associated with 
the undeformed parabola, and i and j associated with the best -fit parabola. 

The radius vector R A from the apex of the original parabola to A' is given as 

= ( X A + 0*1 + ( ax A + V)Jl (3) 

The closest point to A 1 on the best-fit parabola is denoted B (Figure 1) and has the 
coordinates [x#, (a + b)x 2 (3] in the ( X , y) coordinate system. The radius vector from the 
origin of the original parabola to j3 is given as 

Rb = A*1 +Aii + x B i + (a + b)x%] (4) 

where (3\ and are the coordinates of the origin of the best-fit parabola. Denoting the 
angle between the X, and X axes as (3% we have 

*i = * cos /?3 - j sin #3 (5a) 

*2 = sin /? 3 + j cos fo (5b) 

Using Eqs. (5) we can obtain from Eqs. (3) and (4) the error V of A 1 with respect to the 
best-fit parabola 

^ [(ft “h b^)x g V a\j (6) 
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where (xAiVa)’ coordinates of A' in the best-fit coordinate system are given as 


x A = {x A + 0 c ° s (h + (ax\ + j]) sin p 3 - ft cos ft - ft sin ft (7a) 
y A = ~(x A + 0 sin ft + (arr^ + if) cos ft + Pi sin ft - p 2 cos ft 

(7b) 

The point P on the best-fit parabola is found by minimizing \\Rb ~ R'a II res P ec t 
to Xp. In the present work this is done with a Newton-Raphson iteration using X a as an 
initial guess (Xg is the solution to a cubic equation). 

The parameters b, P \ , ft, ft are found by minimizing the root-mean-square (rms) dis- 
tortion 

"rms =2 1 j h V ' VdX ~ ^ ( 8 ) 

where ±h are the limits of the parabola, C; are quadrature weight and X{ are points where 
the deformed parabola coordinates are given. In this work the minimization was performed 
by a conjugate-gradient method using finite-difference derivatives. 

Instead of performing this nonlinear analysis it is standard practice to linearize the 
problem. First we set COS ft = 1, sill ft = ft and neglect higher order terms in Eqs. 
(7) (assuming ft P \ , P 2 , ft are small) to get 


XA = X A + £ + ax\p 3 - p, 
Da = ~x A p3 + ax\ + T) - p 2 


(9a) 

(9b) 
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Next we use a linear approximation to the minimum distance as follows: We set Xq = 
Xa in Eq. (6) and take only the component of V normal to the best-fit parabola at 
Xg — Xa- This normal n is given as 


n = 


j — 2(a + b)xAi 
\Jl + 4(a + 6) 2 x A 2 


( 10 ) 


Neglecting higher order terms the normal component of V can be written as 


v n = v ■ n — —v 0 + fa 


( 11 ) 


where 

Vo = {q - 2ax A f)lJl + 4 c?x\ 


(12) 


and 


f = [x 2 a , —2 ax a, l,xa + 2a z x\]/ \J 1 +4a 2 x\ 


(13) 


The rms error is now defined as 


Vl 



To minimize it we differentiate Eq. (14) with respect to Ot to obtain 


1 

2 h 



v„ 


dVn 

daj 


dx = 0 



Using Eq. (11), Eq. (15) becomes 


Aa = f 


(14) 


(15) 


(16) 
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where 



n 

U T dx = c di( x i)^f {xi) 

i = 1 

n 

v 0 ldx = ^2 CiV 0 (xi)l(xi) 

i=i 


(17a) 

(17b) 


The matrix A is almost singular so that small deformations 7] can result in large 
values of f3\ (the x-translation) and (3$ (the rotation). We can minimize V^ms an 

additional limitation on the size of Of of the form 

a T a < £ (18) 


and this leads to a system of equations 


(A + A I)a = f 


(19) 


where A is a Lagrange multiplier (chosen to satisfy Eq. (18)) and I the identity matrix. 


Best-Fit Paraboloid 

The derivations for the paraboloid parallel the derivations for the parabola given in the 
previous section.. The undeformed shape of the paraboloid is given as 

Z\ = ap\ (20) 

in a coordinate system shown in figure 2. The distortions in the /Cq, 9\ and Z\ directions 
are given by £,T] and respectively, so that point A in Figure 2 moves to position A 1 . 
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The best -fit paraboloid is given as 


(21a) 

(216) 


z = cp 2 
c = a + 6 


The radius vector R ^ from the apex of the original paraboloid to A! is given as 


r 'a = [( Pa + 0 COS 0A-V sin 9 A ]i x 

+\{pA + 0 sin 0A + V cos 6 A ] ji + {ap\ + 0*1 

= x' A h + y' A j\ + z' A ki (22) 


The closest point to A! on the best-fit paraboloid is denoted B (Figure 2) and has the 
cylindrical coordinates [/9g, @b , (a + b)p 2 B ] in the (#, y, Z ) cordinate system. The radius 
vector from the origin of the original paraboloid to B is given as 


Rb — PiH + p2j\ + Pzk\ + Pb cos 9 si + Pb sin 9 b j + (a + b) p 2 B k (23) 


where now Pi , P 2 and ps are the coordinates of the apex of the best-fit paraboloid in 
the , yi, Z\ ) system. The relationship between the unit vectors in the original and 


best-fit systems is given as 



where 



1 

O 

O 

l 


cosp 5 

T = 

0 cosPi —sinPi 


0 


0 sinp 4 cosP 4 


_ —sinPs 



0 

sinp 5 ' 


’ COsPq 

— sinP 6 

O' 

1 

0 


sinP 6 

cosPq 

0 

0 

cos/3 5 _ 


0 

0 

1 

( 
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and /? 4 , (3$ and (3q are rotations around the axes X\, t/i, Z\, respectively. 


Using Eq. (25) we can obtain from Eq. (22), (23) and (24) the error V of A' with respect 
to the best-fit paraboloid. 


v = R b - B! a = ( p B COS e B - x A )i + (p B sin 0 B - y A )j + ( cp% - z A )k ( 26 ) 


where ( X A , y A , Z A }, the coordinates of A' 


in the best-fit coordinate system, are given as 


X A 

v_a 

Za 



( 27 ) 


r ^ It _ MO 

The point B on the best-fit paraboloid is found by minimizing || V | with respect to 
pB and Ob - Doing so we obtain 

0 b = » a +4> ( 28 ) 


where 


tancf) 


y A cos 9 A — x A sin 0 A 


x A cos 9 A + y A sin 9 A 

and p B is the solution of the cubic equation 


( 29 ) 


2 c 2 Pb + pb( 1 - 2 cz A ) - (x A cos 0 B + y A sin 0 B ) = 0 ( 30 ) 


which is closest to The parameters are found by minimizing the root- 

mean-square (rms) distortion 

h 27T ^ 

is = ^2 J J V-VpM - J2 C ^ 2 (p^ e 0 ( 31 ) 

0 0 1=1 
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where h is the limit of p for the paraboloid. As in the case of the parabola the minimization 
was performed by a conjugate-gradient method using finite-difference derivatives. As in 
the case of the parabola we consider also a linear analysis setting COS Pi = 1 , sin/3i = Pi 
for i = 4,5,6. The linear approximation to the minimum distance is obtained by a 
similar procedure to the two-dimensional case: We set ps — Pa an< J @B — @A in Eq. 
(26) and take only the component of V normal to the best-fit paraboloid at pg = p A and 
0 B = $ A . The normal fl is given by 

k — 2 cp A cosO A i — 2cp A sin9 A j 

vA + 4<M 

Neglecting higher-order terms the normal component of V can be written as 

u n = u.n = — Uq + £ l a (33) 

where 

= (C - 2 ap A £)/ yfi + (34) 

and 



f - [ p\,-2ap A cos6 A ,-2ap A sin0 A , 1 , p A sind A (l+2a 2 p\), -p A cosQ A (l+2a 2 p 2 A )\ 


lyJl + ±a 2 (? A 


(35) 


and 


= [b,P\,P2,P3,{34,Ps] 


(36) 
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Note that, as expected because of axial symmetry, 7 ] and do not influence the error. 
The rms error is again defined as 

h 27r 

"rms = ^jj ""VpM ( 37 ) 

0 0 

and the minimization again leads to a system of linear equations. 


Act = f 


(38) 


where 


h 2tt 


A— f f U'pdpdB = 

A A L—\ 


0 0 

h 27 r 


/ 


= J J votpdpdQ = CiVo(pi,6i)£(pi , 0j] 

0 0 L=1 


(39a) 

(396) 


Results for Best-Fit Parabola 

To illustrate the problems associated with finding the vector O' which defines the best-fit 
parabola consider a distortion of the form 

( = (c( 1 - COS (40) 

with 7J — 0. A parabola with a/h = 0.2 corresponding to focal length over diameter ratio of 
0.625 was used. The best-fit parabola was calculated for a very small disturbance £c/h = 
0.005. The best-fit parabola corresponding to this distortion was calculated three different 
ways: 
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(i) By using conjugate gradient optimization procedure to minimize V^ ms based on the 
nonlinear expressions in Eqs. (6) and (8). The resulting surface error is denoted V^i. 

(ii) By solving the linearized problem Eq. (16). The corresponding linear approximation to 
the surface distortion is denoted Vi. 

(iii) By solving the size-limited problem, Eq. (19) for various values of A. 

The results are summarized in Table 1. The first line shows the results obtained with 
conjugate gradient minimization of the nonlinear expression for the error. It is seen that 
the rms value of the error can be reduced by a factor of three. However, there is great 
amplification of the disturbance with the normalized translation (3\ / h being equal to 
0.1571. The linear analysis based on Eq. (16) yielded similar values for the components 
of O'. However, because of the large values of these components the prediction of that 
linear analysis was erroneous. The predicted rms value of 4.9 X 10 -1 compares with a 
nonlinear value of 6.37 X 10 ~ 3 . Thus while the linear analysis predicted a reduction of 
the initial rms by a factor of 3 the nonlinear analysis predicted that the best-fit linear 
parabola actually increased the error by a factor of three and a half. 

The next three lines in Table 1 include size-limited solutions obtained from Eq. (19) 
with various values of A. It is seen that as A is increased the size of a decreases so that 
the linear and nonlinear predictions become close. However, this is accompanied with 
substantial increase in the best-fit rms error. The last line in the table shows a 3- variable 
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solution obtained by setting (5\ to the apex x-translation (zero here). This solution is close 
to the large- A solution from Eq. (19). 

Table 1 shows that we have a dilemma in the construction of a best-fit parabola. Linear 
analysis requires that we eliminate one of the variables or restrict the size of the solution. 
These limitations, however, substantially increase the error rms of the now ‘not-so-best-fit’ 
parabola. The alternative nonlinear solution is complex and costly. 

This type of difficulty is not encountered when the distortion does not require (3\ and 
for its correction. As an example consider a distortion of the form 

£ = & sm — (41) 

The results of the nonlinear and linear solutions are shown in Table 2 for a substantial 
value of £ s /h — 0.04. It is seen that there is hardly any difference between the linear 
and nonlinear solutions. 

Results for Best-Fit Paraboloid. 

Similar results were obtained for the best-fit paraboloid for Cljh — 0.2. For example, 
a distortion of the form 

£ = £ s sin(7r/9/h)sin# (42) 

was considered, and the results for £ s //l = 0.001 are summarized in Table 3. The first 
column shows the results of the nonlinear analysis coupled with the conjugate gradient 
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minimization. The rms error is reduced by about a factor of 7, however there is again 
amplification of the distortion due to the ill-conditioning of the problem with /? 2 /h = 
0.09. The linear analysis shown in the second column produces very similar solution, 
predicting even better reduction in rms (about a factor of 10). However, when the nonlinear 
solution is analyzed using the nonlinear analysis we find that the error actually ijricreaseA 
by a factor of 3. 

The next three columns in Table 3 show the size- limited solutions based on Eq. (19). It 
is seen that, as we put more and more stringent limits on the magnitude of the solution, 
the agreement between the linear and nonlinear solution improves. However, much of the 
reduction in the error is lost, so that we have a ’not-so-best -fit-paraboloid’. Similar results 
are obtained by setting j3\ and /3 2 equal to X and y translation of the apex (zero for the 
example) and solving a reduced 4- variable problem. 

While this dilemma of how to calculate the best-fit-paraboid is difficult, there is a bright 
side to it. The linear analysis gives a reasonable idea of the magnitude of error reduction 
possible with the best-fit paraboloid. 


Concluding Remarks 

An investigation of alternatives for calculating the best-fit paraboloid to a deformed 
paraboloid surface was investigated. In particular we focused on the ill-conditioning as- 
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sociated with the translations perpendicular to the axis of the paraboloid. It was shown 
that this ill-conditioning results in disturbance amplification so that small deformation 
can result in large translations and rotations for the best-fit paraboloid. It was also found 
that eliminating the two translations or restricting their magnitude may result in large 
increases in rms errors. 

The amplification of translations and rotations for the best-fit paraboloid results in 
grossly inaccurate prediction by linear analysis of the rms error. However, the lincai 
analysis may be less inaccurate in predicting the achievable reduction in rms error. 
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Table 1: Best-fit parabola with various fitting schemes, £ c //i = 0.005, initial error 


Vorms/h = 1.75 X 10 3 


Fitting scheme 

b/h 

Pi/h 

th/h 

Pa/h 

u L /h 

vnl / h 






rms 

values 

4- variable 







nonlinear 

.00104 

.1571 

.00426 

.05854 


4.90X 10 -4 

linear 

0 . 

.1574 

0 . 

.05844 

4.90xl0~ 4 

6.37X 10 -3 

\ — < 
o 
o 
o 
o 

II 

0 . 

.1238 

0 . 

.05844 

4.90X 10 -4 

3.90X 10 -3 

A = .0001 

0 . 

.0426 

0 . 

.0146 

8.91 X 10“ 4 

9.74X 10 -4 

A = .0002 

0 . 

.0248 

0 . 

.00785 

9.89X 10 -4 

9.95X 10 -4 

A = .0005 

0 . 

.0113 

0 . 

.00267 

1.07X10 -3 

1.07X10 -3 

3-variable 

0 . 

0 . 

0 . 

-.00163 

1.13X 10 -3 

1.13X10 -3 


Table 2: Best-fit parabola with various fitting schemes, = 0.04, V 0 rms/h = 

8.69 X 10 3 


Fitting scheme 

b/h 

Pi /h 

02 /h 

0a/h 

v L /h 

L / h 






rms 

values 

nonlinear 

.0142 

0 . 

-.00228 

0 . 


6.02X 10 -3 

linear 

.0141 

0 . 

-.00225 

0 . 

5.69X 10 -3 

6.02X 10 -3 
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Table 3: Best-fit parabola with various fitting schemes, ( , s /h — 0.001, initial error 


Vorms/h = 4.868 x 10- 4 


Fitting 

scheme 

nonlinear 

6- variable 

A = 0 

A = 5 x 10- 

linear 

- 6 A = 2 x 10" 5 A = 5 x 10“ 

4- variable 
-5 

b/h 

2.33 xlO -4 

0 . 

0 . 

0 . 

0 . 

0 . 

/?l/h 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

A/ h 

-0.09055 

-0.09055 

-0.06281 

-0.03280 

-0.01686 

0 . 

A A 

1.41X 10 -3 

0 . 

0 . 

0 . 

0 . 

0 . 

A/h 

-0.03366 

-0.03366 

-0.02313 

-0.01174 

0.00569 

7.12X 10~ 4 

A/* 1 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

t'l/h 


5.19X 10 -5 

1.12X10 -4 

2.14X 10 -4 

2.70X 10 -4 

3.29xl0~ 4 


6.91 Xl0 -5 

1.47X 10 -3 

7.05X 10~ 4 

2.79X 10 -4 

2.73X 10 -4 

3.29xl0~ 4 
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